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Abstract 

Penguin contributions, being not negligible in general, can hide the in- 
formation on the CKM angle a coming from the measurement of the time- 
dependent B^{t) TT+vr^ CP-asymmetry. Nevertheless, we show that this 
information can be summarized in a set of simple equations, expressing a as 
a multi- valued function of a single theoretically unknown parameter, which 
conveniently can be chosen as a well-defined ratio of penguin to tree ampli- 
tudes. Using these exact analytic expressions, free of any assumption besides 
the Standard Model, and some reasonable hypotheses to constrain the mod- 
ulus of the penguin amplitude, we derive several new upper bounds on the 
penguin- induced shift \2a — 2aefr|, generalizing the recent result of Grossman 
and Quinn. These bounds depend on the average branching ratios of some 
decays K^K^, K^TT^) particularly sensitive to the penguin. On the 

other hand, with further and less conservative approximations, we show that 
the knowledge of the Ktt^ branching ratio alone gives sufficient infor- 

mation to extract the free parameter without the need of other measurements, 
and without knowing \ Vtd\ or \ Vub\- More generally, knowing the modulus of the 
penguin amplitude with an accuracy of ~ 30% might result in an extraction of 
a competitive with the experimentally more difficult isospin analysis. We also 
show that our framework allows to recover most of the previous approaches in 
a transparent and simple way, and in some cases to improve them. In addition 
we discuss in detail the problem of the various kinds of discrete ambiguities. 

^Laboratoire associe au Centre National de la Recherche Scientifique - URA D00063 
E-mail: charles@qcd.th.u-psud.fr 



1 Introduction 



In a close future, several collaborations — BaBar, BELLE, CDF, CLEO, HERA-B — 
will hopefully make the first measurements of CP-violation in the system [Q . The 
most important consequences concerning the Standard Model (SM) would be the 
determination of the Unitarity Triangle (UT). However, if the measurement of the 
UT angle sin 2(3 seems to be straightforward from both experimental and theoretical 
points of view thanks to the very clean B — > J/'^Ks decay, the extraction of a 
from the standard mode B -k^-k^ is still an open problem |. Since it has been 
pointed out that QCD — and mixed QCD/electroweak — radiative corrections (called 
"penguins") induce potentially large theoretical uncertainties on this angle 0, many 
papers have been devoted to this subject 0. 

In a pioneering paper [Q, Gronau and London have shown that the knowledge 
of the B{B) — s> TT+TT^, 7r°7r°, vr^vr^ branching ratios leads to the determination of 
the gluonic penguin effects, assuming isospin symmetry and neglecting electroweak 
penguin contributions. Then, with this information and the usual mixing- induced 
CP-asymmetry it is possible to get a up to discrete ambiguities. The main drawback 
of this interesting method is the expected smallness of the B — * vr'^Tr" branching 
ratio (10^^-10^^) due to colour-suppression. This fact, combined with the detection 
efficiency of the final state and the needed tagging of the flavour of the i?-meson, 
constitutes a difficult challenge to e^e~ i?-factories and an almost impossible task 
for future hadronic machines — LHCb, BTeV. 

Then it was realized by Silva and Wolfenstein that by extending the flavour 
symmetry to SU(3) one can gain further information on penguin effects, the key point 
being the Ktt modes where the ratio penguin/tree is certainly greater than 1. Con- 
sidering the crudeness of the assumptions made in the original paper in addition to 
SU(3), the method has been extended until a high level of sophistication by several 
authors 0. As a consequence, it is not clear to what extent such complicated geo- 
metrical constructions, plagued by multiple discrete ambiguities, are sensitive to a 
and to the unavoidable theoretical assumptions. Therefore these strategies will give 
conservative results only when a better understanding of non-leptonic 5-decays is 
available. In addition, two simpler SU(3) approaches concerning a have been pro- 
posed by Buras and Fleischer |]^ and Fleischer and Mannel |p respectively, which 
will be discussed in more detail below. 

One can also use a model — usually factorization — to estimate the penguin ampli- 
tude, and then compute the difference between a at the input and at the output. 



as Aleksan et al. H and Ciuchini et al. ITO] did, or directly get a model-dependent a 



as was proposed by Marrocchesi and Paver |[TT| 0. 

Thus, after having hunted [|12|, trapped |jl3| and made the zoology of the pen- 



guin, it is time to begin taming it. To accomplish this task, we first remark that most 



^Throughout this paper, B stands for a Bd meson. 

^ActuaUy we wiU see that the Marrocchesi-Paver method is essentiahy the same as the 
Fleischer-Manuel M one, although the theoretical input is different. 
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of the authors cited above have computed the observables — branching ratio and CP- 
asymmetries — as functions of the theoretical parameters — QCD matrix elements and 
CKM factors, including the angle a. We will follow the opposite way, and show that 
it is indeed a fruitful approach. Although fully equivalent to the "traditional" one, 
it leads to a very important and simple new result: it is possible to express indepen- 
dently of any model 0, and in an exact and simple way all the theoretical parameters, 
including the angle a, as functions of the experimentally accessible observables and of 
only one real theoretical unknown. The latter can be chosen as, e.g., \P/T\, the ratio 
of "penguin" to "tree" amphtudes (which are unambiguously defined below) . It is also 
possible to use as the unknown a pure QCD quantity, free of any dependence with 
respect to \Vtd\ or |V^*^| contrary to the parameter \P/T\; in the latter case, we give 
polynomial equations directly expressed in the (p, rj) plane. We have exploited these 
exact analytic expressions to derive several new and simple results and to recover 
some of the previous approaches. The main points of this paper are: 

• Using the exact parametrization in terms of \P/T\, it is possible to represent 
the information given by the time- dependent CP-asymmetry in the {\P/T\, 2a) 
plane. Of course without any further assumption on the magnitude of \P/T\ 
there is no way to constrain a. But this {\P/T\, 2a) plot provides a nice trans- 
parent presentation of experimental data, where our ignorance of the strong 
interactions is relegated to a single parameter. 

• As soon as one is interested in quantifying the size of the penguin — and indeed 
we are, sin 2a is not a good parameter. One should simply use 2a instead. 
Actually using sin 2a rather than 2a is not wrong, but one loses half of the 
information as we will see in detail below. This is already true at the level of 
the parametrization in terms of \P/T\, and this is also true for all the methods 
allowing to remove the penguin effects, which give generically 2a rather than 
sin 2a, up to discrete ambiguities. To make clear this point which up to now 
has remained confused, we will treat explicitly the example of the Gronau- 
London isospin analysis. On the contrary, the observables depend only on 2q; or 



equivalently on tana, and thus the a ^ vr + a ambiguity is always present |]T4| 



Bounding the magnitude of the penguin allows directly to bound the shift of the 
CKM angle a from the directly observable agg. This can be done using informa- 
tion from decays particularly sensitive to the penguin. For example, assuming 
SU(2) isospin symmetry and neglecting electroweak penguin contributions we 
are able to derive two bounds depending on BR(i? — > vr^vr^) , one of which being 
the Grossman-Quinn bound |T^ while the other is new. Assuming the larger 
SU(3) symmetry, we obtain two new bounds depending on BR(i? K^K'^) 



*In this paper, "model-independent" means "not relying on a particular hadronic model which 
describes non-perturbative physics" . On the contrary, we will assume that the SM holds for the 
parametrization of CP-asymmetries and amplitudes. 
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and A^BR(i? K^ti^) respectively which, not surprisingly, may be more con- 
straining than the SU(2) ones, and which need some, but not all, the usual 
assumptions concerning the neglect of annihilation and/or electroweak penguin 
diagrams. As far as the branching ratios of the penguin-sensitive modes are 
concerned, these bounds do not need flavour tagging and are still valid when 
only an upper limit on the branching ratios is available. In addition, they can be 
slightly modified to be used when the actual value of the direct CP-asymmetry 
in the B —>■ vr+vr" channel is not available, as it is shown below. Depending on 
the actual values of the branching ratios, the theoretical error on a constrained 
by these bounds could be as large as ~ 30° or as small as ~ 10°. In particular, 
the most recent CLEO analyses of the tt+tt^ and K^n'^ modes allow us to 
give for the first time the following numerical bound 

|2a-2aefr| < A, with 25° < A < 59° (1) 

assuming rather weak hypotheses in the SU(3) limit (see § |5.2|) and BR.{B —>■ 
TT+TT") > 0.4 X 10~^ in addition to the CLEO data. 

• Finally, after having stressed that only one hadronic parameter has to be es- 
timated by the theory in order to get a, we give one new explicit example: 
assuming SU(3) and neglecting annihilation and electroweak penguin diagrams, 
we show that BKlB"^ —>■ Ktt^) gives sufficient information to solve a degree- 
four polynomial equation in the (p, 1]) plane, which roots can be represented 
as curves in this plane. Contrary to the Fleischer-Mannel proposal 0, ours 
does not need the knowledge of \Vtd\ or \Vub\, and requires only the measure- 
ment of BR(i?^ — > Kn^) in addition to the usual time-dependent B tt+tt 
time-dependent CP-asymmetry. Alternatively, the knowledge of the modulus 
of the penguin amplitude (or the ratio of penguin to tree) with an uncertainty 
of ~ 30% should provide a rather good estimation of a. This kind of strategy, 
although affected by potentially large theoretical uncertainties, may be neces- 
sary when the more conservative bounds are too weak to be really useful in 
testing the SM. 

The paper is organized as follows: in Section |^, we summarize the main results of 
this work — this section should be of immediate use for the reader not interested by the 
development. In Section ^we fix our notations in writing the general parametrization 
of the amplitudes. With the help of the recent CLEO measurements of non-leptonic 
charmless -B-decays, we give some rough orders of magnitude of the expected penguin 
pollution. Then we derive the equations giving the theoretical parameters, including 
a, as functions of the observables and the theoretical unknown, treated first as a 
free parameter, and latter eventually constrained under reasonable hypotheses. For 
example in Section | we show how to use in our framework the information coming 
from the B — >■ 7r°7r° and B^ vr^vr" decays, to obtain the Grossman- Quinn bound 
and a new similar isospin bound. In Section ^ we exhibit two new bounds, based on 
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the SU(3) assumption, which may be more stringent than the two isospin bounds. 
Then in Section ^ we discuss an exphcit example where the theoretical unknown 
is actually estimated rather than bounded. A reasonable knowledge of a can be 
expected even if one allows a sizeable violation of the theoretical assumptions. In 
Section |^ we discuss how to incorporate and improve some of the previous approaches 
in our language, and clarify some points which have been mistreated in the literature, 
in particular the problem of the discrete ambiguities. Our conclusion is that although 
the penguin-induced error on a is expected to be quite large in the B — > tt+tt^ 
channel, it should be under the control of the theory. Therefore the generalization of 
the methods presented here to other channels is very desirable to get more constraints 
on a. 

This paper has two technical appendices: the first one explains how we got 
the values of the observables from a naive calculation, in order to numerically il- 
lustrate our purpose before experimental data is available and the second one (0), 
following Grossman and Quinn |]15|, shows explicitly the existence of bounds which 
are independent of the measurement of the direct CP-asymmetry. 



2 Summary 

2.1 Exact Model-Independent Results 

Defining the Standard Model B tt+tt^ amplitudes 

A{B^^n+n-) = VudV:,M''^^ + VuV:,M^''^ 
the time-dependent B^{t) tt+tt^ CP-asymmetry 



flcp (t) = Odir COS Am t — Jl — a^jj, sin 2acs sin Am t , 



and the average B, B ^ f, f branching ratio 

1 r 



B 



fj 



BR(i? ^ /) + BR(i? ^ /) 



(2) 
(3) 



(4) 



(5) 



we prove in the following that the Standard Model predicts very simple relations 
between a and \P/T\, \P\, \T\ and 6 = Arg(PT*) respectively, these relations de- 
pending only on the observables 8^^+^^-, Odir and 2aefj and being completely free of 
any assumption on hadronic physics: 



COs(2q; — 2Q;efr) 

\P\' 



1 a^jj. 



1- Jl- a^ij. cos 2aeflf 





p 


2- 


) 


T 





COS 2a 



1 — Jl — a^jj. cos(2a; — 2aefr) 



(6) 
(7) 



4 



1^1 
tan 5 



COS la. 



1 - cos 2aeff 
Odir tan a 



fljij. [cos 2Q;efr + tan a sin 1acs\ 



(8) 
(9) 



Rather than |-P/T|, |P| and |T| which incorporate respectively a iVfd/V^*^!, |Vtd| 
and \y^^ factor, one may prefer to write Eqs. (|-|D in terms of \M^^ IM^'^\ |MW| 
and |M(")| respectively (see the definition (^). As and \y^h\ depend also on the 
UT, it is not possible to express such relations as functions of a alone; instead we 
use the Wolfenstein parametrization and find three polynomial equations in the (p, 77) 
plane. With the definitions of the following combinations of observables 



1 - cos 2aeff , -D^ 



1 — a?: sin 2aefT 



(10) 



and the theoretical parameters IM*^*)! and |M("^| normalized to -\/H~k~/|AV"c 



ch\ 



Rp 
Rt 



\XVr 



cb\ 



R 



T 



WVr 



)|2 



one has two degree-four polynoms depending respectively on Rp/Rx and Rp 



(l-^e)(l-|^)/ + 2(l-Z^.)(l-f^)pV 



(l-^e)(l-f^)^^ 



-2(1 - D,){1 - 2^) - 2D, p\ - 2(1 - D,){1 - 2^) pr/^ - 2D, r]' 
+ (1 - D,)(l - 6^) p2 + 2D,pr/ + [1 + - 2(1 - Z^e)|^] 



+ (l-^c)|^(4p-l) 







'12) 



(1 - D,) p" + 2(1 -D,- Rp) + (l-D,- 2Rp) r]^ 

-2(1 - D,) p3 - 2Ds p^v - 2(1 -D,~ 2Rp) prf - 2D, r/^ 

+ (1 - D,) p2 + 2D, pr^ + (1 + D, - 2Rp) r/^ = , (13) 

and one linear equation depending on R^ (the ± sign being related to a discrete 
ambiguity) 

^1 - D,(p - 1) ± ^2Rt - 1 + D, r/ = 0. (14) 

Eqs. ( p!^ - |I^) are another way of writing Eqs. (|-H) by replacing respectively \P/T\, 
\P\ and |T| by the ratios Rp/Rt, Rp and Rt- the advantage is that the latter 
parameters do not depend on the badly known CKM matrix elements \Vtd\ and \ V*ij\. 
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2.2 Phenomenological Applications 

It has become standard in the CP-hterature to use several phenomenological assump- 
tions, some of which can be very good while some others can be strongly violated. 
As a result, it is often not easy for the reader to know exactly which approximations 
are used by the authors, and thus to make his own opinion about the accuracy of 
these theoretical prejudices. In this paper, we will try to state clearly what kind of 
hypotheses we use in addition to the SM; some of the results that we derive rely on 
a few reasonable assumptions chosen in the list below. 

• Assumption 1 \P/T\ < 1. This very conservative bound should be distin- 
guished from the small penguin expansion. 

• Assumption 2 SU(2) isospin symmetry of the strong interactions. 

• Assumption 3 SU(3) flavour symmetry of the strong interactions. 

• Assumption 4 Neglect of the OZI-suppressed annihilation penguin diagrams 
(cf. Fig. I). 

• Assumption 5 Neglect of the electroweak penguin contributions. 

• Assumption 6 Neglect of the VusV*^^ contributions to the —>■ K^n^ ampli- 
tude. 



Upper Bounds. We have found several quantities bounding the shift of the true 
2a from the experimentally accessible 2aes, among which ([T6| ) is the Grossman- Quinn 
bound [|T5|1, while the others are new: 



if sin 2acfT >0 0<2a<27r — 2 arcsin(sin 2acs) 

if sin2aofr < — 2 arcsin(sin 2aefr) < 2a < 2it [assuming 1], (15) 



2a — 2aes\ < arccos 



1 - al. 



dir 



1 - 2 



[assuming 2 and 5], (16) 



2a — 2aes\ < arccos 



1 — a 



dir 



1-4- 



[assuming 2 and 5], (17) 



2a — 2apff| < arccos 



1 a^jj. 



[assuming 3 and 5], (18) 
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\2a — 2Q;efr| < arccos 



1 - 2W 



2 ^K^-n"^ 



[assuming 3 and 4]. (19) 



Following Grossman and Quinn |jT5|| we show that, under the same hypotheses, 
the above upper bounds still hold if one replaces in Eqs. ( p!5| - |T9|) adir by zero and 2aes 
by 2aes where the latter effective angle is defined by 



sign(cos 2Q;efr) = sign(cos2a;) , sin2Q;ofr 



For example, one has the bound 



2a — 2aes\ < arccos ( 1 — 2A 



B-JT+Ti 



1 — a^jj, sin 2aes 



[assuming 3 and 4] 



(20) 



(21) 



and so on. Although these bounds on \2a — 2aea\ are weaker than the ones on 
1 2a — 2aefT|, they have the advantage that they do not depend on the measurement 
of adir^ indeed, the angle 2aeff is accessible (up to a twofold ambiguity) from the 
sin Am t term only (see Eqs. (^) and (^)). Therefore the experimental uncertainty 
should be smaller for 2aes than for 2acs- 

In addition to these upper bounds, we derive a lower bound on |2a — 2aefr| in 
^ 5.41 to which we refer the reader for more details. 



Determination of a. We propose a new method for the extraction of a — up to 
discrete ambiguities, which improves the Fleischer-Mannel proposal [^. 

The idea (often used in the literature) is to estimate the modulus of the penguin 
contribution with the help of the Kir^ decay. We avoid the problem of 

knowing \Vtd\ by using directly the polynomial equation ([T3|) in the (p, 77) plane, with 
the theoretical parameter Rp given by 

JO 

Rp = X^-JSlL [assuming 3, 4, 5 and 6]. (22) 

This typically leads to draw four allowed curves in the (p, t]) plane, which in the limit 
i?P — > 0, reduce to the two circles representing the no-penguin solution sin 2a = 
sin 2aefT. 



3 Theoretical Framework 

3.1 Standard Model Parametrization of the Amplitudes 

The aim of this paragraph is to recall some already known results and to fix the 
notation used in this paper. 
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The time-dependent rate for an oscillating state B^{t) which has been tagged as a 
meson at time t = is given by (for simplicity the e~^^ and constant phase space 
factors are omitted below |) 



TT TT 



1 cos Am t — Im —AA sm Am t . 

2 2 \P 



(23) 



r {B\t) 

where 

A = A{B^ ^ TT+TT-) , A = A{W TT+TT-) (24) 
and q/p = exp(— 2z/5) in the Wolfenstein phase convention, which provides an expan- 



sion of the CKM matrix in powers of A = \Vus\ ~ 0.22 With this convention, 
one has 

/5 = Arg(-l-;,) , 7 = Arg(-\C,) , Arg(y,,) = ^(A^) , (25) 

while the other CKM matrix elements are real (up to highly-suppressed A" terms) 
and the angle a is given bya = 7r — /? — 7. Defining 



Odir 



1 r 



Ul 



BR(5° 

Ul 



TT TT 



BR{B^ 



TT TT 



) 



2«efT 

the rate (ESf) becomes 



|A|2+ 

Arg (-AA* 



(26) 
(27) 

(28) 



r (B%t) TT+TT - 



1 -I- fldir COS Am t — a/ 1 — a^jj. sin 2aefr sin Am t 



(29) 



The time- dependent CP-asymmetry reads 



acp(t) 



r (SO(t) ^ TT + TT-) - r (fiO(t) ^ TT 



TT 



r(50(t) ^TT+TT-) + r(50(t) 



TT^TT 



= adir COS Am t — yl — a^^. sin 2aes sin Am t . 
We may define another effective angle by | 



sign(cos 2Q;efr) = sign(cos 2a) , sin 2acs = — (Mir sin 2a, 



such as 



c^cp (t) = Odir COS Am t — siu 2aes sin Am t . 



(30) 



(31) 



(32) 



^Tiny differences between phase space of the various channels discussed in this paper are 
neglected. 

^As the sign of cos2aeff is not observable, it can be defined arbitrarily. However, the exact 
definition is important for the derivation of the bounds (cf. Appendix . 
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adir is the direct CP-asymmetry, while sin 2aes ~ y ^ ^ '^dir sin 2aofr is the mixing- 
induced CP-asymmetry. In the absence of penguins, one has adir = and sin 2aes = 
sin2aefr = sin 2a. The experiment allows the measurement of three fully model- 
independent observables, one is CP-invariant {8.,,+^,-), while the two others are CP- 
asymmetries (ttdir and sin2aefr, or adir and sin2aefr)- 

It has often been assumed in the literature that the dominant penguin amplitude 
is the top-mediated one, with the consequence that this amplitude is proportional 
to VtdV^l. This assumption has received a lot of attention recently |]18], |1^. In any 
case, the w-penguin contributions which are proportional to V^dl^*6 (as well as other 
contributions such as exchange diagrams) can be incorporated in the definition of the 
"tree" amplitude. Contributions proportional to V^^V^, — "charming penguins" [|T^ — 
can be rewritten, by CKM unitarity (VcdV*)^ = —VudV*^, — VtdV^^), in terms of the two 
other combinations. Thus, just from the weak phase structure of the SM, we may 
write the — > tt+tt" physical amplitude as 

A = VudV:,M^''^ + VtdV,lM^'^ 

= e+*^T^+^- + e-*^P^+,- (33) 

and similarly for the other 27i channels [| 

^^0^0 = A(fiO ^ ttV) = -L(e+*^T,o,o + e-''^P^o^o), (34) 

= A(5+ ^ 7r+7r°) = -^(e+*^T^+,o + e-*^P^+^o). (35) 

Let us stress that there is absolutely no approximation in writing Eq. ( p3H35| ): Tj^t, and 
P„-^ are CP-conserving complex quantities, defined by the weak phase that they carry, 
and they incorporate all possible SM topologies such as trees, penguins, electroweak 
penguins... In this sense, many (not all, however) of the methods proposed previously 
for the extraction of a in the top-dominance assumption are in fact still valid if the 
latter hypothesis is relaxed f\ The CP-conjugate channels are obtained by reversing 
the sign of the weak phases: 

As said above, neglecting penguin diagrams (P^^r = 0) gives 



2a,s = Arg ^T,+^^T;+^_e-2^^ ] = 2a . (37) 




From now on we will denote in the whole paper 

Ttt + TT' = T , Pn + TT- = P , (38) 

^Note that P^+^o comes from electroweak penguins, and/or from isospin symmetry breaking. 
®0f course, numerical estimates of quantities like \P/T\ may be greatly modified by charming 
penguins pi, pi. 
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and by economy of language, we will call the T and P amplitudes "tree" and "penguin" 
respectively, although T gets contributions from u- and c-penguin diagrams ^ 

In this paper, we will also consider the B — > K'^K^, B — > K^tt"^ and B^ — > Kt:^ 
decays. For the former, we adopt the following parametrization 



A{B^^K^K^) = 



(39) 



while for the latter it is convenient to expand on the CKM basis (KtsKTfe) KsKfe) (recall 
that VcsV*i^ is real in the Wolfenstein convention): 



A{B^ 
A{B+ 



K+7T' 



K + TT- 



TT 



(40) 
(41) 



As far as the B — >■ K^K^ amplitude is concerned, we have used the notation Tj^j^ to 
make apparent the resemblance with the other channels; however it should be stressed 
that this decay is a pure penguin process. Actually Tj^j^ represents the contribution 
of the long-distance u- and c-penguins |^ . 

Let us repeat that Eqs. (p3|-p6|) and ( p9| - p| ) rely only on the Standard Model. 



3.2 General Bounds 

we will denote by Bf j the CP-conserving average branching 
1 



Similarly to Eq.( 
ratio 



For example. 



B 



fj 



BR{B ^ /) + BR{B ^ /) 



(42) 



Bktt^ — - 



BR{B^ 



BR{B- 



TT 



(43) 



and so on. 

From the discussion in § |3.1| , it is clear that the SM predicts each B^ ^ f decay 
amplitude as the sum of two terms carrying two different CP- violating phases 0i, 02 : 



A{B^ ^ f) = e+'^'Mi + e 



+i<t>2 



2 , 



A(50 ^ f)=r]f (^e-'^'Mi + e-^'^Wa) (44) 



where Mi and M2, although complex numbers, are CP-conserving and the sign r]f 
depends on the CP of the final state, e.g. rjf^n^n^) = +. Thus the average branching 
ratio (B2f) writes 



B 



fJ 



|Mip + |M2p + 2Re(MiM*)cos(0i -02) 



(45) 



^Be careful that our definition of "tree" and "penguin" amplitudes, relying on CP-phases, is 
slightly different from the one used in Refs. although the consequence is the same: these 

so-defined amplitudes are unambiguous and physical quantities. 
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Note that we express the amphtudes squared in "units of two-body branching ratio" . 
For fixed values of M2 and 4>i — 4>2, ^jj as a function of Mi takes its minimal 
value when Mi = — M2cos(0i — 02); in which case i3jj|mm = |M2psin^(0i — 02)- 
Interverting the roles of Mi and M2, we obtain the following (exact) general bounds: 

|Mi|2sin2(0i -02) < Bfj, |M2psin2(0i -02) < Bjj. (46) 

Such inequalities have been previously employed for the demonstration of the so- 
called Fleischer-Mannel bound . We will use Eq. (|46| ) extensively throughout this 
paper. 



3.3 Orders of Magnitude 



The recently updated CLEO analyses of i?-decays into two light pseudoscalars give 
precious information on the various quantities discussed in this paper [IC 



< 0.84 X 10"' [90% CL] , 
= (1.4±0.3±0.2) X 10-^ 
= (1.4±0.5±0.2) X 10"^ 



(47) 



Thanks to this experimental information, it is possible to derive a crude lower bound 
for the ratio \P/T\. Indeed, from Eqs. (Bl) and (E6|) we have 



ITI sin a < By, 



(4J 



Furthermore, while the tttt penguin is proportional to VtdV^l, the Kn penguin is 
proportional to VcsV*fj (cf. Eqs. (0) and (^0| - pD ). Thus we have 



IPI 



td 



cb 



X 



(49) 



where Eg is a typical scale of the B Ktt branching ratios, assuming that these 
channels are dominated by the QCD penguin, and that the QCD-part of the penguin 
matrix elements are of the same order for vrvr and Kn. Thus we obtain 



> A| sin7| 



(50) 



where we have used |Vtd/Kf,| = A sin 7/ sin a. 

Numerically, from the current SM constraints on the (p, //) parameters [^, we 
have I sin7| > 0.6. The CLEO data (0) suggest BJB^+^~ > 1/0.84 which gives 



> 0.14. 



(51) 



Note that contrary to some claims, factorization gives typically \P/T\ ~ 0.15 [§] and 
is not ruled out by the CLEO data yet, although it is only marginally compatible. 
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Although this calculation is only illustrative, it is clear that the penguin contri- 
butions pose a serious problem for the extraction of a from the CP-asymmetry. More 



complete analyses show that the ratio |-P/T| can easily be ~ 30% or even 40% [19 
Therefore one can not avoid to define theoretical procedures allowing to reduce the 
penguin uncertainty, or at least to control it. This is the subject of the present paper. 

Let us add here one comment on the size of the direct CP-asymmetry. With a 
ratio |-P/T| of order (9(20-30%), it is straightforward to show that a strong phase 
Arg(PT*) of order (9(10-20°) is sufficient to generate a direct CP-asymmetry of 
order (9(10%) (cf. Appendix |^). While perturbative calculations of such phases 
predict in general very small values it is likely that non-perturbative effects will 
considerably enhance these Final State Interactions (FSI) Q [^. Thus we expect 
that sin2aefr and Odir will be measured with comparable statistical accuracy |p6 |. 

At this stage, we refer the reader to Appendix 0, where we calculate the relevant 
parameters and observables in a naive way, in order to numerically illustrate the 
phenomenological results that we derive below. 



3.4 Some Exact One-Parameter Results 

Let us first consider only the B^ii) — > tt+tt^ rate, Eq. (|29|) . As pointed out above, 
there are three observables: the average branching ratio and the CP-asymmetries 

i3^+7r- , o-dir and sin la^si ■ (52) 

From sin2aeff, one gets 2aeff up to a twofold discrete ambiguity. While the vanishing 
of the penguin amplitude P implies la^si = 2a, the SM description of the rate (pQ]) 
involves four parameters, three of which are CP-conserving while the fourth is the 
CP- violating angle a (see Eqs. ( |33D and (^)), namely 

|T|, |P|, 5 = Arg(PT*), and a, (53) 

one overall phase being irrelevant, and after the use of q/p = exp{—2if3). Thus the 
presence of the penguin amplitude forbids the measurement of a, the number of 
parameters being greater than the number of observables. However, as we will see, 
it is possible to express a in terms of the three observables and of one of the four 
parameters. The latter can be chosen as either \P/T\, \P\, |T|, or 6. 
From (|33D and (^), we deduce 

-{2isma)P = e-'^A-e'^A, (54) 
(2isina)T = e'^A-e-'^A, (55) 

"'^'^Note also that in the — > oo hmit, smce u- and c-penguin can contribute, such phases between 
T and P amphtudcs arc in principle 0(1). Indeed as a perturbative calculation suggests the 
real and imaginary parts of the long-distance penguins are of the same order. 
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that can be rewritten as 



{2i sin a)P 
{2i sin q;)T 



P 



a-1a 

P 



(56) 
(57) 



Calculating the ratio |-P/Tp from ( ^^ - |57| ) and using the definitions ( ^^^ - |28| ) we obtain 
the very important, although very simple, equation: 



cos(2a — 2a;efr) 



1 



1 a^jj. 



1 



1 - a\-^ cos 2acfr 





p 


2 - 


) 


T 





(5J 



Thus Eq. ( [J^ j defines 2a as a four-valued function of\P/T\: indeed, 2aes is known 
up to a twofold discrete ambiguity, while Eq. ( ^8|) also contains a twofold discrete 
ambiguity as long as 2a is concerned, because of the cosine function. Note that these 
two discrete ambiguities are of a different nature: the 2aes — > vr — 2aes ambiguity is 
inherent to CP-eigenstate analyses, while the 2a — 2aes —{2a — 2aes) ambiguity 
is generated by the penguin contributions. One can also view the latter ambiguity 
by saying that the no-penguin solution 2a = 2Q!eff appears as a double root of the 
general cosine equation (0), which degeneracy is lifted by the penguin contributions. 
Let us add here one comment on the discrete ambiguity generated by getting 



28 



q) and 



a from 2a, that is the a n + a ambiguity. From 
one sees that the three observables Odir and 2acs are invariant under the 

transformation a^7T + a,6-^TT + 6. Thus these observables depend only on 2a 
and 26, or equivalently on tana and tan 5. Without any further assumption on the 
strong phase, the a — > tt + a ambiguity is irreducible , the signs of sin a and sin S 
being related by the equation: 



sign(sina) x sign(sin5) = sign(adii 



(59) 



As far as the SM is accepted, the latter ambiguity not a real problem because the 
constraints on the UT select only one of the two ambiguous solutions — we already 
know that < a < vr because rj is positive |^ — and moreover because these solutions 
cannot merge as they are separated by tt. This is not the case, obviously, for the 
ambiguity 2acs — * vr — 2aeg. In the following we will always express our results in 
terms of 2a and 26, or tana and tan 5. 

It is also possible to derive very simple relations expressing the parameters \P\, 
\T\ and tan 5 as functions of \P/T\, or equivalently, as functions of 2a which is itself a 
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function of \P/T\ through Eq. (^81). Indeed, from (|26| - p8|) , (|53|) and (^6]-| 



I I ^ 

tan 5 



1 — cos 2a 

1 — COS 2a 



1 - a/ 1 - a^i^ cos(2a - 2aefr) 



1 — a^jj, cos 2q;, 



eflf 



Odir tan a 



1 — Jl — a^jj, [cos 2aes + tan a sin 2aeff ] 



we get 0: 
(60) 
(61) 
(62) 



Note the following important point: Eq. ( |60D (resp. Eq. (|6lD) gives 2a as a four- 
valued function of \P\ (resp. |T|), and Eq. (|6^) gives 2a as a bi-valued function of 
25; this feature puts the parameters \P/T\, \P\, \T\ and 2S on an equal footing: each 
of these is candidate to be the single theoretical input. 

Let us discuss the no-penguin limit of Eqs. ( [58| ) and ( |60| - p^ ): \P\ ^0 (and 



0). In this limit, Eqs. 



and 



reduce simply to 2a = 2aefr, as 



thus Odir 

it should. Eq. ( pT| ) reduce to |Tp = i3^+,r- independently of a, as also expected. 
Finally Eq. (|62D becomes indefinite in this limit, because 6 itself becomes indefinite. 
Note the important point that the parameters \P/T\ and \P\ measure directly the 
size of the penguin and thus of the shift \2a — 2aes\ while the parameters |T| and 6 
carry only poor information on the size of the penguin (for example the no-penguin 
relation |Tp = i^Tr+vr- can still be verified with a non- vanishing \P\, for particular 
values of the parameters). 

Lastly, we stress that Eqs. (58-62|) are fully equivalent to the original Eqs. (|33D 
and ( ^6]) together with the definitions (|26| - p8D and (^31). They are exact relations 
between the theoretical parameters taken two by two, depending only on the ob- 
servables. The most obvious application of this result is the modelization of the 
theoretical error induced by penguin effects on the CKM angle. For example the 
reader may choose his favourite model or his favourite assumptions to estimate \P/T\ 
as well as its associated error. This range of values of \P/T\ propagates into four 
cleanly-defined, although model-dependent, discrete solutions for 2a and their theo- 
retical errors, thanks to Eq. ([581) . The same can be done using, instead \P/T\ and 
Eq. (IID, the parameters \P\ and Eq. or \T\ and Eq. (|61|), or 26 and Eq. (^) []. 
As the model is used to calculate only one real quantity, this procedure should be safer 
than the ones proposed in, e.g., Refs. and |Ty|. We will give practical examples 
of this strategy in Sections || and ^. 

Another application, which is more conservative but less informative, is to use 
channels where the penguin may be dominant and, when related by a flavour sym- 



^^The fact that only the weak angle a is present in Eqs. (|58| ) and (|6(]|-|62[) originates from the 
peculiar SM prediction that the i?" — mixing is dominated by the top-loop, just cancelling the 
CP-phase of the penguin defined by Eq. (^3|), and is not related to the dominance (or not) of the 
top in penguin loops. This is obviously not the case, e.g., for the decay B Kstt^ where both a 
and (3 enter in the game. 

^^Note that the extraction of a from Eq. (|6^) would require a very accurate and unlikely knowledge 
of ,5. 
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metry to the parameter |P|, can help to bound the shift \2a — 2a;efr|. This is the case 
for B — >■ 7r°7r°, K^K^ and K'^n'^, as explained in Sections ^ and ^. But before that, 
we shall insist now on the model- independent features of Eqs. (|58D and ( |60| - |6^ . 



3.5 Plotting the CP- Asymmetry in the ( ^ , 2aj Plane 



From I cos(2a 
obtained: 



2«cfr)| < 1 and from (|58D the following allowed interval for |P/Tp is 



1 



''dir 



1 



< 



a^jj. cos 2aes 



< 



1 + 



1 a^jj. 



1 — a% cos 2aeff 



(63) 



The lower bound on |P/Tp is induced by the direct CP-asymmetry — it becomes triv- 
ial in the limit Odir 0; indeed, if the latter is non zero, then a non-vanishing penguin 
amplitude follows. As the sign of cos 2aes is not observable, (p3| ) defines actually two 
different intervals for \P/T\, one for the branch corresponding to cos2aefr > and 
the other for the branch corresponding to cos2aefr < 0. 

Assuming Odir and sin2aefr have been measured, Eq. (^) allow to plot 2a as a 
function of \P/T\ varying in the interval (|63|): on Fig. [l| we represent the two distinct 
branches corresponding to the two possible signs for cos2aefr- 



cos 2(Xeff < 




COS 2(Xeff > 



2 . 5 



3 

\P/T\ 



-100- 



Figure 1: The CKM angle 2a as a function of\P/T\ obtained from Eq. ([5^, fora^h = 
0.12 an(i sin 2Q;eff = 0.58 (see Appendix^ for a description of this numerical example). 
The solid curve corresponds to cos2aefr < and the dashed one to cos2aeff > 0. Note 
that \P/T\ varies in the range ^6^ ) which depends on the sign of cos2aes. For a given 
value of \P/T\, there are in general four solutions for 2a. 



Furthermore, the three remaining equations ( |60| - |6^ ) together with 2a given by 
Fig. allow to represent |P|/v^HJh^, \T\/^yB^r+^T- , and 25 as functions of \P/T\ 
(Fig. I). 

Let us summarize the main properties and virtues of Eqs. (^) and (|60| - |6^ ) and 
of Figs. |I^: 
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0.5 1 1.5 2 2.5 3 

\P/T\ 

150v 

100-- 
2S (deg.) : 
50v 



0.5 1 1.5 2 2.5 3 

\P/T\ 




1V5- 2 2i5 ■;;,::;3„„.. 

\P/T\ 



Figure 2: The \P\, \T\ amplitudes and the strong phase 25 respectively as functions 
of \P/T\ obtained from Eqs. ( |5^ j and ^6di - \6^ ), for the numerical example a^iT = 0.12 
and sin2aeff = 0.58 (see Appendix^). The solid curves correspond to cos2a;efr < 
and the dashed ones to cos2aefr > 0. Note that \P\ and \T\ diverge when 2a 0, 
\P/T\ —> 1 because the CP- asymmetries a^n and sin2a;efr are kept fixed: \P\ and \T\ 
have to be very large to produce CP-violation with a small CP-phase a. 



• These are absolutely exact results, relying only on the SM, and provide a nice 
representation of what kind of model-independent information can be obtained 
from the measurement of the time-dependent B^{t) — ^ 7r"'"7r~ CP-asymmetry 
only. The {\P/T\,2a) plot may be used to present the experimental results 
which hopefully will be available from i?-factories in the next years. 

• Fig. ID shows that the linear approximation in \P/T\, which is used in some 
papers |] is indeed very good for \P/T\ < 1 as far as 2a — 2aes is concerned. 
However, Eqs. (|58|) and ( |60l - ^ ) expanded to first order in \P/T\ give expressions 
which are not particularly simpler, and thus it is more convenient to keep the 
exact formulae. 

• Eqs. (^) and (|5U|-|5^) are not invariant under the transformation a ^ — a. 
This does not properly mean that the a ^ — a ambiguity is lifted: the sou- 
venir of this ambiguity lives in the invariance of Eqs. (|5^) and (|60| - |62|) under 
des ^ f ~ "^eff because the sign of cos 2aes is not known. This means however 
that sin 2a is not a good parameter, indeed the penguin effect is not the same for 
the solutions corresponding to cos 2aes > than for the others corresponding to 
cos2aefr < 0, as Fig. |l| clearly shows. In particular, the solutions corresponding 
to cos 2acs < are more affected by the penguin uncertainty which is an impor- 
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tant information 0. As long as the penguin is not strong enough to change the 
sign of cos2aefr/ cos 2a, we actually expect cos2aefr < from the current SM 
constraints on the UT [23]. To be illustrative, let us plot sin 2a as a function 
of \P/T\ using Eq. and compare with 2a as a function of \P/T\ (Fig. 0). 
Obviously, four curves in the (|P/T|, sin 2a) plane gives half as less information 
than four curves in the (|P/T|,2a) plane. If we reconstruct the \P/T\ 2a 
curves from the \P/T\ — > sin 2a ones, we will get eight curves among which 
four are "wrong" solutions (Fig. We discuss further this point in § 7J^, with 
the explicit example of the Gronau-London construction. We conclude that 
one should not express the penguin effect in terms o/ sin 2a — sin2aofr as it is 
sometimes done in the literature [^, 




(a) 



(b) 



Figure 3: (a) sin 2a as a function of \P/T\ obtained from Eq. for the numer- 

ical example Odir = 0.12 and sin2aefr = 0.58 (see Appendix^). The two branches 
corresponding to the two possible signs for cos 2a cs o,re not degenerate, (b) 2a as a 
function of \P/T\ obtained by computing arcsin(sin2a) and vr — arcsin(sin2a).' the 
comparison with Fig. || shows that the dashed curves are wrong solutions. 



• Bounding the absolute magnitude of the penguin directly allows to bound the 
shift 1 2a — 2aofr| (and vice-versa) thanks to Eq. ( |5^ ) or Eq. (pO|). For example, 
the very conservative estimate \P/T\ < 1 (Assumption 1) leads to the simple 
bound 

cos(2a — 2aefr) > cos2aefr . (64) 

Of course, this bound does not allow a precise measurement of a. Nevertheless, 
with only a very weak assumption, it provides an allowed interval for 2a (we 
have taken into account that the sign of cos2aefT is not known): 

if sin 2aefr >0 0<2a<27r — 2 arcsin(sin 2aefr) 

if sin2aefr < — 2 arcsin(sin 2aefr) < 2a < 27r [assuming 1]. (65) 

As explained in Appendix ^, this bound implies a weaker one, obtained by 
replacing above 2aeff by 2aeff where the latter effective angle is defined by 

'This has already been noticed by Gronau in Ref. [ p7[ . 
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Eq. (^). The only advantage in using 2aes instead of 2aes is that the former 
angle follows directly (up to a twofold ambiguity) from the sin Ami term in 
the time-dependent CP-asymmetry ( p2D independently of Odir; thus the experi- 
mental uncertainty on 2aes is expected to be smaller than on 2aes- 

For the numerical example that we have chosen (Appendix ^), sin2aefr = 0.58, 
we obtain from the bound (^) < 2a < 290°, which is of course not very 
informative. 



3.6 Exact One-Parameter Polynoms in the {p, 77) Plane 

For evaluation by hadronic models, or by using phenomenological assumptions such as 
in Section |], the theoretical parameters \P/T\, \P\ and |T| may be not suitable as they 
depend on QCD matrix elements times |Vtd/Vu*b|, \Vtd\ and \V*i,\ respectively. Indeed 
the latter CP-conserving CKM factors are badly known, therefore they would intro- 
duce an additional uncertainty in combination with the theoretical model-induced 
error for the estimation of the QCD-part of the matrix elements. In the literature, 
this problem has been solved by scanning the whole allowed domain for (p, rj) , 
by simply assuming that such CKM factors would be known from other measure- 
ments [^, or by expressing jV^d/KTbl as a function of a and P, the latter angle being 
determined from future CP-measurements in the B — > J/'^Ks channel ITTI . 



We think however that it is more convenient and more transparent to decouple the 
different and intricated problems related to the determination of the UT. Fortunately, 
such an attitude is simple to handle with, thanks to the CKM mechanism which 
predicts strong relations between CP- violating and CP-conserving quantities: indeed 
the SM says that a, \Vtd\ and \V*i\ are functions of (p, r/) [p!7| : 

a = Arg ( j-P-''l \ ^ (66) 

\Vu\ = AlKfel x\l-p-zr]\, (67) 

Kbl = MVcbl X \p + tv\. (68) 

The above relations, inserted in Eqs. (|58|) , (l60| ) and (|6l|) , permit us to reexpress the 
latter as equations in the (p, t]) variables, depending on the theoretical parameters 
|M(*) |/|M*^") I, \M^^^\ and |M*^")| respectively. Indeed we define the following combi- 
nations of observables 



= V 1 - fldir cos 2aeff , D s = \J 1 - a^,. sm 2aes , (69) 
and we introduce \M^^^\ and IM*^"-*!, normalized to v^B^r+^/l-^Kfel (recall their def- 



inition (p3D, and that the amplitudes squared are in "units of two-body branching 
ratio" ) 



Rp_ 




2 


Rp = 






5 



|MW|2 _ \M(^ 



)|2 



'cb\ 



R^=\XV,,\'^ (70) 
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For simplicity, we neglect the uncertainty on Vcb, and take Vud = Vtb = 1- 
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Then Eqs. (|58D and (|60|) are respectively equivalent to the following degree-four poly- 
nomial equations, the first depending on Rp/Rt only and the second on _Rp only: 




(1 - D,) + 2(1 -D,- Rp) p'rf + {l-D,- 2Rp) r]^ 

-2(1 - D,) p3 _ 2D, p^T] - 2(1 -D,- 2Rp) pr/^ - 2D, r]^ 

+ (1 - D,) p2 + 2Ds p7] + {l + D,- 2Rp) rf = Q. (72) 

When M*^*) = (the no-penguin case: Rp = and thus ajir = 0), Eqs. ( |7ID 
and ff73) reduce to 



[1 — cos 2aeff] 



P 



1 



V 



sin 2aofr 



1 



1 2 



2(1 — cos 2a, 



eS) 



2(1 — cos 2a, 



off J 



(73) 



which is the equation squared of a circle. This circle is just, as expected, the one 
defined by 2a = ATg[—{l — p — iri)/{p+iri)] = 2acs and can be obtained geometrically, 
by using the definition of the UT and solving the equation 2a = C^^. Actually, 
the sign of cos2aefr is not known and we get in fact two circles. When M*^*) 7^ 0, 
each of these two circles splits into two curves — this splitting is reminiscent of the 
2a — 2acs —{2a — 2acs) ambiguity of Eqs. (|58D and (|60D : the no-penguin case — 
the circle — appears as a double root of the general case — a degree-four polynomial 
equation, as already noticed above when discussing Eq. (|5^). 

Likewise Eq. (|6T|) is equivalent to the following linear equation, depending on Rjn 
only: 

'1 - (p - 1) ± J2Rt -1 + D,r] = 0. (74) 



The ± sign is reminiscent of the 2a —2a ambiguity of Eq. (|6]J 

As the parameter Rt does not know much about the size of the penguin, the 

no-penguin limit of Eq. ([7^) is not particularly interesting. 

The important feature of Eqs. (|7T]-|73|) is that the parameters Rp/Rx, Rp and 

Rt — defined by Eqs. ( ^3]) and (ffOD — are pure QCD quantities times |AVcfep/i3^+7r-, 

i.e. they can be expressed as matrix elements of the Weak Effective Hamiltonian 

times known factors. 

Thus the reader may choose a pure hadronic model to estimate Rp/Rx, Rp or 

Rt, and report it in Eqs. (0), ( |72D or ([7^ ) respectively, then getting a polynomial 
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equation which roots, represented as curves, summarize the domain in the (p, 1]) 
plane which is allowed by the measurement of the time-dependent i? — > vrvr CP- 
asymmetry. Some examples of this strategy are given in Section ^ where we use 
some phenomeno logical assumptions to estimate Rp, and in Section |^, where we 
suggest to improve the proposals of Fleischer and Mannel P] and Marrocchesi and 
Paver |jTl| by solving the problem directly in the (p, 1]) plane. 



4 Using Isospin Related Decays 

In this section, we will assume: 

• SU(2) isospin symmetry of the strong interactions (Assumption 2). It is well 
known that this flavour symmetry is indeed very good; in any case, the violation 
of SU(2) should be completely negligible compared to the theoretical errors 
discussed in this paper. 

As the Effective Weak Hamiltonian is a linear combination of AI = 1/2 and AI = 3/2 
operators, one has the triangular relations (remember the notation (|33| - ^ ) and (^8|), 
and see |,|T1): 

Ttt+ttO = T + TttO^O , -Ptt+ttO = P + PttO-it" ■ (75) 

As the QCD-penguins are pure AI =1/2 amplitudes, the Pn+wO amplitude come only 
from electroweak penguins. Thus we define Pew = Pn+wO to get 

A{B^^n^n-) = e^^T + e'^^P, (76) 



1 



A{B' ^ ttV") = [e^-^ + e-^^(PEw - P)\ , (77) 
A(P+ 7r+7r°) = i= [e'^ (T + T^o^o) + e'^^PEw] • (78) 
The second assumption we will make here is: 

• Neglect of the electroweak penguin contributions in B tt+tt^, tt^tt^, tt^tt^ (As- 
sumption 5). That is Pew = in Eqs. ([76| - |78D . The problem with this approx- 
imation arrives when considering B — > where, on naive grounds (short 
distance coefficients and factorization of the matrix elements), the electroweak 
penguin — which is here colour-allowed — is not particularly negligible. However 
the repercussion on the extraction of a is expected to be negligible E^, B^] , and 



in any case smaller than the gluonic penguin effects. See also the discussion in 



In the framework of these two assumptions, Gronau and London have shown that 
the knowledge of the B{B) tt+tt^, vr^vr^, vr^vr^ branching ratios in addition to the 
time-dependent B^{t) tt+tt^ CP-asymmetry leads to the clean extraction of a, up 
to discrete ambiguities. In § \! .1[ we reexpress the Gronau-London isospin analysis in 
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our language. In particular, we clarify the problem of the discrete ambiguities, which 
up to now has remained confused in the literature. 

Unfortunately, it is well known that the isospin study might be experimentally 
difficult to carry out, if the B tt^tt^ mode is as rare as expected because of colour- 
suppression. Therefore alternative methods have to be developed. 



4.1 Two Upper Bounds on \2a — 2aes\ from B 



In Ref. Grossman and Quinn have derived an interesting bound on the shift 
|2a — 2acfT| (Eq- (2.12) of their paper) which takes in our notation the simple following 
form: 

\2a - 2acs\ < arccos (l - 2-^) . (79) 



This bound derives from the isospin relations ([f^-|78|) and from the geometry of the 
Gronau-London triangle (see 0, |13[) when the electroweak penguin amplitude, i.e. 
Pew, is neglected. The physical meaning of this bound is simple: the B tt^tt^ 
branching ratio cannot vanish exactly unless both the tree and the penguin amplitudes 
in vr'^vr^ vanish, in which case 2a = 2aeff in tt+tt". 
As explained in Ref. 
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this bound is useful when the B n n rate is too 
low, in which case only the average branching ratio Bt^o^o (that can be obtained from 
untagged events only), or even only an upper bound on this quantity, is available. 
Thus, either the B — »• tt'^tt^ channel is strong enough to allow a full isospin analysis, 
or the rate is indeed very small and bounds the penguin- induced error on a. 

It is not difficult to derive the bound (|7^) in an analytical way, different from the 
geometrical approach of Ref. ; here we give only the main line of the demonstra- 
tion. Neglecting the electroweak penguin (Pew = 0) in Eqs. ( [76l - [7^ ) we can form 
the ratio 8^,0^,0 / Bt,±t^o and consider it as a function of the complex parameter T^+ttO- 
Minimizing this ratio with respect to the latter parameter gives the inequality: 



LP sm a < — X 



B^ 



Then 2a is constrained thanks to Eq. 



1 2a — 2acs\ < arccos 



''dir 



1 - 2 



Bt,±^0 



[assuming 2 and 5]. (81) 



As < Odir < 1, the above bound is slightly more stringent than the bound (|79|) , and 
reduce to the latter when adir = 0. 

Under the same isospin symmetry and neglect of electroweak penguin hypothe- 
ses, it is straightforward to derive another similar bound, not given in the original 
paper ||15|. Indeed, using the general bounds (|46|) for the penguin in Eq. ( ^4]) gives 
simply: 



iPMsin^a < 2B^o^o 



iS2) 
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where the factor two is related to a Clebsh-Gordan coefficient, i.e. to the wave 
function of the vr*^ meson and the Bose symmetry. Once again we use Eq. (BDI) to get: 



\2a — 2aea\ < arccos 



1 '^14 



1 - al^ V ^-+-- 



[assuming 2 and 5]. (83) 



To which extent the Grossman-Quinn bound ( pTD is better than the bound (pSf ) 
or vice-versa depends on the actual values of the branching ratios 2B.„±t,o vs. i^Tr+Tr-: 
in fact, neglecting penguin and colour-suppressed contributions would lead to the 
equality 2i3^±^o = ^B^+tt-, while the factorization assumption, predicting a construc- 
tive interference between the colour- allowed and colour-suppressed contributions in 
the —>■ 7r^7r° channel, tends to favour (^) compared to ( |83| ) 0. A technical 
advantage of the bound (|83|) over (|8T1) is that it does not require the measurement 
of ;B^±^o, which may be less well measured than 3^^+^^- because of the necessary 7r° 
detection, and because it is expected that Bt,±t,o < Bt,+t,-. Note in passing that the 
Grossman-Quinn bound follows from the isospin constraints on both tree and penguin 
amplitudes, while our bound comes only from the isospin constraints on the penguin. 
Of course it can be checked that the two bounds are fully compatible, in the sense 
that when the bound ( pi]) is saturated then the bound (|83|) is automatically satisfied 
and vice-versa. 

As shown by Grossman and Quinn |[T^, and as redemonstrated for consistency in 



Appendix g the above bounds imply weaker ones, obtained by replacing Odir by zero 
and 2aeff by 2aefr where 2aes is defined by Eq. (0). Thus we have: 

\2a - 2aes\ < arccos (l - 2^^] , (84) 

\2a — 2acs\ < arccos 1 — 4— I [assuming 2 and 5J. (85) 

As already stressed, the advantage in using 2aes instead of 2aes is that the former 
angle follows directly (up to a twofold ambiguity) from the sin Am t term in the time- 
dependent CP-asymmetry (|32D independently of Odir, and thus does not require the 
measurement of the latter. 



5 Using SU(3) Related Decays 

In this section we will assume a larger flavour symmetry, namely 

• SU(3) flavour symmetry of the strong interactions (Assumption 3). One could 
argue that such an assumption should not be too bad in energetic two-body de- 
cays, although we know that a typical SU(3) breaking quantity is Ifx — fnl/fn ~ 

"'^^Grossman and Quinn give another bound depending on both S^o^o/S^±^o and 6^07^0/6^+^- 
(eq. (2.15) of their paper |jl^). As it is more compHcated and presumably numerically similar, we 
do not report it here. 
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23%. Actually, our present knowledge does not permit a reliable quantitative 
estimate of such a symmetry breaking in B decays, especially for the penguin 
amplitudes that we are interested in. In any case, our understanding of this 
problem is expected to improve with both theoretical and experimental pro- 
gresses. 



5.1 An Upper Bound on \2a — 2aeff| from B K'^K^ 

Very similarly to the isospin analysis and the B tt^tt^ case, it is possible to de- 
rive a bound on |2q; — 2Q;efr| depending on BK{B K^K'^). Indeed Buras and 
Fleischer ^ have proposed a SU(3) analysis which rely on the measurement of the 
time-dependent CP-asymmetry in the B K^K^ channel to disentangle the penguin 
effects m. B -k^-k^ (see § [7.3| ). However the B —>■ K^K'^ channel is a pure h ^ d 
penguin process, and its rate is presumably rather small (10~''-10~^). Nevertheless, 



due to the bounds (|46|) , -B^oxo cannot vanish unless both T^o^o ^"^^ ^roI^o vanish in 
Eq. (pOf). Thus either Bj^o'xo is large enough to do the Buras-Fleischer analysis, or it 
is vanishingly small and one expects that \2a — 2aefr| is constrained by Bj^oj^ thanks 
to the SU(3) symmetry. Similarly to the case of tt^tt", an upper bound on Bj^o-j^ is 
sufficient to get constraints on a. 

In addition to the SU(3) flavour symmetry introduced above, we need the following 
assumption: 

• Neglect of the electron] eak penguin contributions in B ^ tt^ti^ , K^K^ (Assump- 
tion 5 ). In Ref. 0, Buras and Fleischer argue that this approximation is better 
than its equivalent in the Gronau-London construction. Indeed, in the latter 
case the electroweak penguin is colour-allowed, while in the present case it is 
colour-suppressed. However one has to remind that FSI effects may invalidate 



the notion of colour- suppression |3^ . 
Within the above assumptions, we have 

1^1 = \Pko1^\ ■ (86) 

Then we repeat the demonstration given above for the B ti^ti^ channel to obtain 
from Eqs. (||), (|4|), and (13): 



\2a — 2aes\ < arccos 



1 - a^ir V ^^+^- . 



[assuming 3 and 5]. (87) 



Likewise (Appendix under the same hypotheses there is a bound independent of 
adir, using the angle 2aeS- 

\2a — 2aes\ < arccos ( 1 — 2 J^°^" [ [assuming 3 and 5]. (88) 
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Hence, analogously to the isospin analysis and the bounds and (p3D, our 
bounds ( P7| - |55D may be useful when the B K^K'^ channel is too rare to achieve 
the full Buras-Fleischer analysis, and thus only the value of Bj^oj^, or even an upper 
limit on this branching ratio, is available. 



5.2 An Upper Bound on \2a — 2aes\ from B 



It has been known for a long time that the B — > Kit decays can help the extraction 
of a from the time- dependent B ^ titi CP-asymmetry by constraining the pen- 
guin amplitudes 0. Indeed, the latter are doubly- Cabibbo-enhanced by the ratio 
l^sKfe/l^sKTb)! "with respect to the tree in these Kn decays. However, in addition 
to the unavoidable SU(3) assumption, people are often lead to neglect annihilation 
diagrams and/or electroweak penguins and/or u, c-penguins and/or Final State In- 
teraction (FSI) in expressing the B ^ titt amplitudes in terms of the B Kn 
ones 0. Such ill-defined approximations have been questioned in the recent liter- 
ature 



30, 31 



in connection with the so-called Fleischer-Mannel bound on sin^ 



7 |22 



Here however, in addition to SU(3), we will only use the following approximation in 
comparing Eqs. (^3|) and (^): 



Neglect of the OZ I- suppressed annihilation penguin diagrams (Assumption 4). 
The topology of these diagrams is represented in Fig. ^. We will need to neglect 
these diagrams only for the P-amplitude, i.e. when the quark in the loop is a 
t or a c (recall Eq. (|33|)). When the flavour in the loop is t, the suppression is 
perturbative, due to a linear combination of short- distance Wilson coefficients 
(mfe); on the contrary, the same diagram with a c quark is non- 



which is 



a: 



10 1 . In addition these diagrams are 



perturbatively suppressed by the OZI-rule 
usually expected to be suppressed by the annihilation topology. Thus they are 
probably very small and negligible compared to the SU(3)-induced theoretical 
error. 





(a) 



(b) 



Figure 4: (a) OZI-suppressed annihilation penguin diagram. This diagram is OZI- 
suppressed and is neglected within Assumption 4 because it does not contribute to 
B K^-n^ . (b) Non-OZI-suppressed annihilation penguin diagram. This diagram 
is not OZI-suppressed, although it has annihilation topology. It is not neglected within 
Assumption 4 because it contributes to both B tt+tt^ and B 
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In particular, we do not neglect the electroweak penguin as it produces the same 
amplitude in S — > K'^ir^ and B — >■ tt+tt", assuming SU(3) and neglecting OZI- 
suppressed penguins. Then we get simply from Eqs. (^), (^) and (|40|) 



IPI 



X IP, 



sm a 



^9) 



where the geometry of the UT has been used in writing |VJd/V^{,|. From Eqs. (^OD 

(90) 



and (^), we get the following bound on \Pk+tt-\- 

|P/^+^-|^sin^7 < BK±n^ 
) and dH) we obtain □ 



Combining it with Eqs. 



1 2a — 2acs\ < arccos 



1 - 2X' 



E 



[assuming 3 and 4]. (91) 



Thus this bound is a quantitative realization of the well-known fact that the penguin 
in 7171 is A-suppressed with respect to the penguin in Kti: if X^BKin^/BTr+n- is not 
too large, this means that the penguin cannot be too large in i? — vrvr. Note that 
similarly to the B — >■ tt^tt" and B K^K'^ channels, we have not assumed penguin 
dominance in B ^ Kn although of course the bound should be more interesting 
when the penguin amplitude really dominates in the latter decay. 

From the experimental point of view, the bound that we have found in Eq. ( pT]) 
should be considerably less affected by statistical uncertainties than the bounds (^TJ), 
( p3D and (0). Indeed, rather than measuring the branching ratio of very rare decays 
such as i? ^ 7r°7r° or i? ^ K'^K^^ the use of the bound ( pT]) needs to know Bk±-kT 
which is (9(10^^). For this reason the CLEO data ( ^T] ) already help to give a very 
interesting and non trivial estimation of the r.h.s. of (|91|). Indeed ( pi]) imply \2a — 
2a;eff| < arccos(l — 2X'^Bk±tt^ /B^^+t,-) and from (|47| ) we have 0.81 x 10"^ < Bk±-kT < 



2.0 X 10-5 and B^+^- < 0.84 x 10"^ at 90% C.L.. Assuming furthermore 0.4 x 10"^ < 
Stt+tt- — otherwise the study of CP-violation in the vrvr channel would be very difficult, 
independently of penguins — we obtain the bound 



2a-2apff|<A, with 25° < A < 59° . 



(92) 



Thus, although these data indicate that the extraction of a will not be an easy task, 
they are still compatible with a relatively small penguin-induced theoretical error. We 
would like to stress also that to our knowledge, this is the first time that a numerical 
upper bound on the theoretical error on a is given rather model-independently, with 
only mild theoretical assumptions and before the experimental value of is available 
by itself. It is expected that experiment will give an accurate determination of the 
r.h.s. of ( PID quite soon. Unless we are unlucky and -B^+,r- is much smaller than 

) and 
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This is a somewhat miraculous feature of the SM: sin 7 cancels between Eqs. 
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expected, the theoretical errol on a constrained by the bound (0) should not excess 
~ 30° while it can be as small as ~ 10°. In comparison, the current knowledge of a 
is roughly 40° < a < 140° []. 

Finally one has again a bound independent of adir where loifS is involved: 

1 2a — 2acff| < arccos {\ — 2A^ ''^~~~~\ [assuming 3 and 4]. (93) 

Let us note in passing that the inequality (0), together with the assumption 
l-Pft'+TT-P = BR(_B='= — i> K-K^^ = Bk-k± (that we will use in Section |^), leads to the 
Fleischer- Mannel bound on sin^7 |22|. However, the bound (pOD is exact indepen- 
dently of the assumption \Pk+-k~\'^ = Bk-k±- 



5.3 Upper Bounds on \2a — 2Q;efF|' Numerical Examples 

In Appendix ^ we define a typical set of parameters for the quantities involved 
in the channels that we are interested in. This set of parameters, compatible with 
the CLEO data (^7]) , allows to compute the various observables (branching ratios 
and CP-asymmetries), and in particular to estimate numerically the bounds we have 
derived until now: 



1 2a — 2aefj| < arccos 



1 



1 a^jj. 



33.4° 



(94) 



2a — 2apff| < arccos 



2a — 2apff| < arccos 



2a — 2apff| < arccos 



1 a^ir 



1 a^ii. 



1-4 



1-2- 



1 



1 - 2W 



B 



40.3° 



27.1° 



29.9° 



(95) 
(96) 
(97) 



The true value being 2a — 2acfr = +26.7° (Appendix ^), the bound (0) is very close 
to be saturated for this set of parameters. Note also that the bounds ( p4[]9^) are 
numerically close one from one another, that just follows from our set of parameters 
and needs not to be true in general. As said above, it may happen in practice that the 
experiment gives only an upper bound on the suppressed channels B.„o.„o and Bj^qj^, 
in which case the bound ( pTf ) will certainly be more informative as Bk±-kt is already 
measured and hopefully the ratio Bk±t^t /B^^+t^- will be known with high accuracy 
very soon. 



^^We stress that although there are presently very weak constraints on sin 2a 
case for a itself. 



this is not the 
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In any case and for the illustrative purpose, we will examine the case of the 
bound \2a — 2Q;efr| < 30°, which is in the ball-park of (pi|-P7|). On Fig. ^ we show the 
constraints of such a bound in the {\P/T\,2a) plane, and in the (p, r/) plane. The 
latter are obtained by plotting the circle defined by 2a = C*, which equation is (cf. 
also Eq. (0)): 




(a) (b) 



Figure 5: (a) The bound \2a — 2aes\ < A in the {\P/T\,2a) plane, obtained as 
explained in the text, for the numerical example Odir = 0.12, sin2Q;eff = 0.58 (see 
Appendix^ and A = 30° (solid curves limited by the small dashes), (b) The same 
bound in the (p, rj) plane, obtained by plotting the circle 1 ^^ ) for 2a varying in the 
interval [2aes — A, 2acs + A]. There are two families of circles corresponding to the 
two possible signs for cos2Q;eff . In the background is shown a crude representation of 
the early- 1998 allowed domain 



3i]. 



Thus the bounds (|9J-^^ should give important and rather safe information on 
the angle 2a as it is apparent on Fig. ^ Even if the penguin-induced error on a 
may be large, it is bounded by theoretical arguments, which is already an important 
statement in view of the possible tests of the consistency of the SM. 



5.4 A Lower Bound on \2a — 2aeff | from B K^tt^ 

It is clear from the examples in § |5.3| that a lower bound on |2a — 2aes\ would be a 
valuable information: it would permit to eliminate some region around 2a = 2aes and 
to get four separate intervals for 2a (cf. Fig |I]) instead of the two big ones represented 
on Fig. ^. Thus one may look for a lower bound on the absolute magnitude of 
the penguin. However, without any further theoretical assumptions (see Section |), 
such a lower bound cannot be obtained using branching ratios only (for example the 
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bound discussed in § ^]3| is not theoretically justified). Conversely, using direct CP- 
asymmetry in the -B^ — >■ K^tt"^ decay, it is possible to get a lower bound on |P|, as 
well as a slightly improved upper bound with respect to the bound (|9lD. The idea 
is the following: if a direct CP-asymmetry in the B — > i^'^vr^ channel is detected, 
then it proves that this mode is fed by both tree and penguin contributions. As the 
latter is related by SU(3) to the penguin in the B —>■ tt+tt" channel, and thus to the 
penguin-induced shift \2a — 2aeff|, one gets a lower bound on this quantity 0. 
Analogously to the derivation of Eq. (|60D , we get from Eq. (^) 



2 IP, 



K + T, 



|2 • 2 

I sm 7 



l-Kl'cosC 



B 



where 



7^ 
^dir 



^ BR(P° ^ K+7V-) - BR(P0 _^ K- 



7r 



(99) 



(100) 



BR(PO ^ K+n-) + BR(PO ^ i^-7r+) 

should be relatively easy to measure for this self-tagging mode, and C is a useful 
short-hand for the phase 



C = 27 - Arg [A{B° K+7i-)A*{B^ ^ R-n 
Then, the inequality | cos^l < 1 together with Eqs. ( p9D and (p9| ) imply 



(101) 



Kit] 
dir J 



1- [a 



and from Eq. (|60|) 



2a — 2aes\ < arccos 



l3K±n^ < 2|P|^sin^a < X^l + \/l-[a 



B 



1 ttjij. 



1 - 1 



(102) 



(103) 



arccos 



1 - «dir 



l-AMl-^/l-k^-l'^ 



''dir 



B^ 



<|2a-2aeff|, (104) 
[assuming 3 and 4]. 



Note that if 



A- ( 1 _ ^/l _ \a^<^f ) < ( 1 - ,/l - al^ ] B^ 



dir ; *-^7r+7r 



(105) 



the SU(3) lower bound in ( |102| ) is useless as it is automatically verified thanks to 
Eq. (|58|) and the (exact) bound (^). Thus this lower bound is only useful in the 

"'^^Note that a non-vanishing direct CP-asymmctry in the tt+tt^ ehannel gives aheady a lower 
bound on the penguin contributions through Eq. (^) . However the saturation of the latter bound 



imply only 2a = 2acs. 
stronger than Eq. (p3). 



Thus one should look for a lower bound on the penguin that has to be 
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configuration where the direct CP-asymmetry is very small in the B —>■ tt^tt" channel 
(odir — > 0) but large in the B — > K'^n^ one (it becomes trivial in the limit —>■ 0), 
in which case the inequality ( |105| ) is not verified. As an example, it can be checked 
that the set of parameters defined in Appendix ^ verifies ( |105| ) . However, keeping the 
same branching ratios and choosing the parameters such as a^ir = and = 0.5, 
the bound (|104|) is not trivial: 

8° < |2a-2aefr| , (106) 

while the bound ( |103| ) represents only a tiny improvement over (|9l|) . 

Actually, one easily obtains similar lower bounds from the two previously studied 
channels, namely B tt^tt", B — > K'^K^. However, the experimental detection of 
direct CP-violation in these suppressed channels may be a difficult task. Should it 
be feasible, one may do the full Gronau-London and/or Buras-Fleischer analyses (see 
Section |^). 



6 Using the B"^ Kn^ Decay to Determine Rp 
with Further Assumptions 

In this section, in addition to the hypothesis made in § |5.2| (Assumptions 3 and 4), 
we will assume more specifically that the two following approximations hold (to an 



accuracy to be determined) in Eqs. (^0|-p! 



Isospin symmetry and neglect of electroweak penguin contributions in B ^ 
K^n"^, Kir^ (Assumptions 2 and 5). Note that the isospin symmetry is a 
consequence of the already assumed larger SU(3) symmetry. Neglecting the 
electroweak penguin, which is here colour-suppressed |2^, we are allowed to 



write Pr+tt- = —Pk°-k+ p3|- 



Neglect of the VusV*^^ contribution to the B^ i^°7r+ amplitude (Assump- 
tion 6). That is, T^o^+ = 0. Using a diagrammatic decomposition of the ampli- 
tude , we have T^o,+ = |KsKl,l(^a + M,-Mt) and P^o^+ = \V,sV:j,\{M,- Mt), 
where Ma is the tree annihilation amplitude and M„ (resp. resp. M^) 
is the u- (resp. c- resp. t-) penguin. Note that T^o^+ is suppressed by 
iKtsKTbl/l^csKll ~ 2 X 10^^ compared to the dominant amplitude Pko-k+- Thus 
we have presumably iKsKTbll^^ - < iKsKbll^c - and it is often 
assumed that annihilation processes are negligible due to form-factor suppres- 
sion P, which then lead to |Tft:o^+| <^ |Pr-o^+|. 



It is clear that Assumption 6 is on weaker grounds than the others made until 
now 0. Accepting it nevertheless, one is lead to many applications [§, H among which 



19 



In particular, it implies a non-trivial relation between FSI phases 
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the most recent one is the Fleischer-Mannel bound 



• 2 / ^K±7rT 

sm 7 < 



B 



(107) 



The latter has been recently questioned The problem is that FSI effects may 
invalidate the notion of colour-suppression for the electroweak penguin, thus leading 
to Pr+tt- 7^ —PkOtt+ [0- Furthermore, the same effects may enhance annihilation 
diagrams, involving a significant KtsKT^ contribution to K^tc^ and a possibly 

measurable direct CP-asymmetry in this channel [0. We will not discuss this subject 
here. Rather we stress as previous authors that the B KK decays may help in 
constraining the FSI effects [^, Q . In particular, the very easy to detect B K^K^ 
mode is fed only by annihilations diagrams. CLEO has already given an interesting 
bound on its branching ratio [T^ : 



Bk+k- < 0.24 X 10" 



C.L.l. 



Thus, either the FSI effects are non-negligible and the K^K~ final state should be 
detected very soon, or they are eventually out of reach of experiment and a stringent 
bound on Bk+k- should be obtained As claimed by the authors of Refs. |31[], FSI 
effects may easily invalidate the bound ( |107| ); indeed, to get a significant constraint 
on 7, we need the ratio Bk±-kt /Bk-k± to be sufficiently less than 1 in order to be 
not too much affected by a reasonable theoretical uncertainty induced by the neglect 
of electroweak penguin and annihilation contributions. On the contrary, for the case 
we are interested in, namely the extraction of a, we do not need Bk+ttI' /Bk-k± < 
1, and we will see that even in the presence of a sizeable violation of the above 
assumptions, we can get interesting information in the (p, rf) plane. In other words 
our method concerning a is useful whatever the values of the branching ratios are. 
However the Fleischer-Mannel bound is not affected by SU(3) breaking, while our 
method is. Note also that Fleischer pOf and Gronau have proposed very recently 
extensive methods which may help to control FSI and electroweak penguin effects for 
the extraction of 7. 

Returning to the problem of a, we use the above hypotheses to write IPr+tt-I = 
\Pko^+ \ = \A{B+ K°7i+)\ = \A{B- K^7r-) \ and thus (recall the notations 



\XVr 



cb\ 



B-rr+n- 



> Bk-k^ 



[assuming 3, 4, 5 and 6]. 



(109) 



The above determination of Rp can be used to insert in Eq. (|7^). Of course, the 
reader who does not agree with the assumptions leading to Eq. (|109|) can use its own 
model to estimate Rp. Thus the method described here is very general, and is in any 



^"Note that the most recent CLEO analyses [|T6[ give Bk±itT /BK-n^ ~ 1; thus the bound (107) 
becomes useless, even neglecting the theoretical uncertainties associated with it. 
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case weakly model-dependent as it depends on only one estimated parameter. The 
results shown below in the (p, 77) plane are quite typical of what can be obtained with 
such a method. 

However, at this stage there is still a problem in using Eq. ( [1091 ): it is clear that 



we have to give a theoretical error associated with the above determination of the 
penguin amplitude. As a guess, we will simply allow a relative violation of Eq. (|109|) 
of the order of 30% and 60% respectively (at the amplitude level), and leave for the 
future any justification of these values. Actually, as long as this error is less than 
100%, the method described here is more powerful than the hounds derived in the 
previous sections 0. 

On Fig. ^ we solve Eq. (|7^) with the theoretical input ( |109| ), and with the nu- 



merical values obtained in Appendix ^ Note that with our set of parameters, the 
Fleischer-Mannel bound becomes trivial (sin^ 7 < 1) but it does not prevent to get 
useful results from Eq. ([721). Fig. ^ shows that with a reasonable 30% relative vio- 



lation of the theoretical assumptions (at the amplitude level) leading to Eq. (|109[ ), 



the time-dependent B rnr CP-asymmetry defines a small allowed domain in the 
(p, rj) plane, much more informative than the more conservative bounds derived in the 
previous sections. This statement is quite general: if there is a way to estimate the 
parameter \P/T\ (or \P\ or Rp/Rx or Rp) with an uncertainty of order ~ 30%, then 
Eq. ( pHD (or Eq. (^) or Eq. ( [7TD or Eq. (|T^)) will give rather strong constraints on a 
(or on the allowed domain in the (p, 77) plane). We will see in § and that the 
isospin analysis is not much more better in this respect because it is plagued by more 
discrete ambiguities. Finally we stress that from the experimental point of view, our 
proposal is very favourable: in addition to the usual time-dependent B^{t) tt+tt^ 
CP-asymmetry, our analysis require only the measurement of the Kir^ average 

branching ratio, which is already measured (cf. Eq. (0)). In this sense our proposal 
represents an improvement with respect to the Fleischer-Mannel's M, because the 



latter needs the further knowledge of Bt,±t,o and \Vtd\ (cf. § [7.4| ) 



7 Recovering and Improving some of the Previous 
Approaches 

In this section, we will explain how to recover in our language the Gronau-London , 
the Buras-Fleischer |0, the Fleischer-Mannel and the Marrocchesi-Paver [jlT| pro- 
posals, and in some places we will propose improvements of these methods. 



Unfortunately, it is not clear if the relation (109) is good at less than 100% relative error. 
Model-dependent criticisms do not predict such a huge violat ion of Assumption 6 ||3l|| , however in 
our case we have to take into account SU(3) breaking in Eq. ( |l09|) . 
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0-1 




0-1 



Figure 6: The solutions of the degree-four polynomial equation in the (p, rj) plane 
for the numerical example Odir = 0.12, sin 2aeff = 0.58 and Rp = \^BKn±/B^+^- = 
0.061 (see Appendix ^). A guess value for the relative theoretical uncertainty on 
Rp in Eq. HlOdj ) is assumed: respectively 30% (a) and 60% (b), at the amplitude 
level. There are four families of curves corresponding to the two possible signs for 
cos2aefr, and to the cosine discrete ambiguity of Eq. ( [g^j which is hidden in the 
polynom (fT^j. In the background is shown a crude representation of the early-1998 



allowed domain j^/ . 



7.1 The Gronau-London Isospin Analysis 

Gronau and London have proposed a clean method to get rid of the penguin-induced 
shift on a |13| by measuring all the B ^ tttt branching ratios in addition to the 
time-dependent CP-asymmetry (^O]). Rather than repeating the geometrical demon- 
stration contained in the original paper, we give here the equivalent analytical for- 
mulae and show the isospin construction in the {\P/T\,2a) plane. 

The Gronau-London method relies on the isospin symmetry of the strong inter- 
actions: after having defined 



$ = ATg{A^+^oA* 
$ = Arg(X-,oA* 

simple trigonometry in Eqs. (|76| - |78D gives: 



cos $ 



cos $ 



1 



V2\A\\A^+^o\ 
1 

V2\A\\A^-^o\ 



;iio) 
111) 



(112) 
(113) 



Eqs. ( |112[]113|) are not yet sufficient to trap the penguin. However, setting Pew = 
in (^) implies Arg (^^A^-^^oAI+^q^ = 2a and thus 

2a = 2acfr + - • (114) 

To summarize, measuring the B ^ txtx branching ratios allows to extract the angles 
$ and $ (up to a fourfold discrete ambiguity which corresponds to the four possible 
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orientations of the Gronau-London triangle |l^) thanks to Eqs. (|112fjll3|) . As the 
CP-asymmetry gives 2aefr up to a twofold discrete ambiguity, it is possible to get 2a 
and \P/T\ from Eqs. ( |114| ) and (0) up to an eightfold discrete ambiguity, as Fig. |^ 
shows. 
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200 + 

150 
2a (deg.) 

lOOf 
50 




Figure 7: The eight solutions of the Gronau-London isospin analysis, for the numer- 
ical example Odir = 0.12, sin2aefr = 0.58, i3^±^o/i37r+7r- = 0.71, 8^,0^0 /Bt^+t^- = 0.061 
and a direct CP-asymmetry in the tt^tt^ channel equal to ajjj,'^ = 0.32 (see Ap- 
pendix^), (a) In the {\P/T\,2a) plane, the dots represent the central values obtained 
from Eqs. ^112^\11^ ), while the solid curves (limited by the small dashes) represent 
the allowed domain when assuming that 2a is affected by a 4° uncertainty due to 
electroweak penguin contributions, as explained in the text, (b) The same allowed do- 



main is represented in the (p, rj) plane, where it is obtained by plotting the circle 
for 2a varying in the eight solution intervals. In the background is shown a crude 



representation of the early-1998 allowed domain [S^j. 



Let us show explicitly that expressing the problem in terms of sin 2a is somewhat 
misleading: from Fig. ^, we can plot the eight solutions of the isospin analysis in 
the (|P/T|, sin2a;) plane, as showed in Fig. ||(a). Now if we forget Fig. ^ and try to 
get the solutions in 2a from Fig. ||(a), we obtain the sixteen solutions of Fig. j^b), 
among which eight are obviously wrong. Note two important points which have been 
mistreated in the original paper |^ and to our knowledge in the subsequent literature: 
First, there are eight solutions in terms of 2a and also in terms of sin 2a 0. Second, 
the isospin analysis determines 2a rather than sin 2a. 



7.2 Defining the Error Due to the Electroweak Penguin 

One may wonder on the size of the electroweak penguin, which is neglected in the 
isospin analysis. Several authors have estimated this contribution, which turns out 

^^If, in addition, the mixing-induced CP-asymmetry in the tt'^tt*' is measured, there are still two 
solutions for 2a (and thus four for a in [0, 27r]), contrary to what is said in Refs. In any case, 

the measurement of this asymmetry is expected to be very difficult. 
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Figure 8: (a) The eight solutions of the isospin analysis in the {\P/T\,sm2a) plane, 
for the same observables as in Fig. j^. (b) The solutions in the {\P/T\,2a) plane 
obtained by computing arcsin(sin2Q;) and vr — arcsin(sin2a).- the comparison with 
Fig. shows that the crosses are wrong solutions. 



to be a few percents of the dominant B — >■ Tr~^7i~ amplitude |2^. If this esti- 
mation is correct, then the corresponding uncertainty on a may be a few degrees, 
which is neghgible compared to the most optimistic simulations of the statistical un- 
certainty In any case, a simple parametrization of the electroweak penguin 
effects can be obtained: indeed, when Pew 7^ 0, there are two new parameters, namely 
Arg(PEW^*) and |Pew|; and one new observable which is the direct CP-asymmetry 
in the P^ 



7r^7r° channel 



BR(P^ 



BR(P- 



TT TT 



BR(P^ 



0-] 



+ BR(P- 



(115) 



which vanishes when Pew ~^ 0. Similarly to the case of the strong penguin, as 
discussed at length in this paper, it is possible to express a as a simple function of 
the observables of the Gronau-London isospin analysis, the direct CP-asymmetry in 
the B"^ vr^TT^ channel and the unknown parameter |Pew|- The same technique 
leading to Eq. ( ^8|) allows to find 



cos(2a — 2agg) 



1 



1 - 



1 



1 



^dir 



COS 2a' 



eflf 



EW 



(116) 



where 2agg is the value of 2a when Pew = (see Eq. (|114|) ) 

2a'^ff = 2aeff + $ - $ . 



(117) 



Thus Eqs. (|116| - |117D describe the departure from the isospin analysis (|114|) due 



to the electroweak penguin contributions. As a particular case, we obtain the bound 



1 2a — 2a'^fi\ < arccos 
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~ 2 















118) 



^^However one should keep in mind that the effect of the electroweak penguin on the B 
branching ratio is not negligible in general. 
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which was derived in Ref. ||28[. Note that the ratio I-Pew/^tt+ttoI is rather independent 
of the size of the colour-suppression, although the impact of |Pew| on BR(i? vr^vr^) 
is not negligible. 

Thus, whatever the way to estimate the parameter IPew/^tt+tt" I? one is lead to a 
simple and weakly model-dependent definition of the theoretical error on a induced by 
the electroweak penguin. Using factorization for the estimation of the r.h.s. of (|118D , 
we find typically \2a — 2a'^g\ < 4°; this error has been reported on Fig. |^ for illus- 
tration. This figure shows the sensitivity of the isospin analysis with respect to the 
discrete ambiguities: with an error on a as small as 2° (here this error comes from the 
electroweak penguin contributions, but unfortunately there are also the uncertainties 
of experimental origin which we have not considered), the four separate solutions 
(for a given sign of cos2aefr) tends to merge quite quickly. This is not a posteriori 
surprising: indeed these four solutions are separated because of the QCD penguin 
contributions, i.e. because of a relatively small effect; they become degenerate in the 
no-penguin limit. Comparison with Figs. ^ and |^ suggests actually that the more 
simpler approaches to control the penguin effects described in the previous sections 
may be competitive with the more complete isospin analysis, unless the observables 
of the latter are known with a very high accuracy. 

The main drawback of the Gronau-London analysis is the expected rarity of the 
B 71^71^ channel, which branching ratio is expected to be about 10^^-10^^. The 
neutral pions are not easy to detect, and one needs to tag the flavour of the S-meson 
in order to get separately |A^o^o| and |y4^o^o|, according to Eqs. ( |112| - |113| ). The 
small number of effectively useful events expected at an e~^e~ 5-factory constitutes a 
difficult challenge to the experimentalists while the impossibility to detect two neutral 
pions in future hadronic machines does not improve the situation. This shows the 
interest of the bounds (ISTD and (E^). 



7.3 The Buras-Fleischer Proposal 

Considering the experimental difficulties associated with the Gronau-London analysis, 
Buras and Fleischer have proposed an alternative way to get rid of the penguin 
uncertainty, using SU(3) and the time-dependent CP-asymmetry of the pure penguin 
mode B K^K^ 0. They argue that the SU(3) breaking effects are of the same 
order as the electroweak penguin uncertainty of the isospin analysis. 

The idea is simple: similarly to the B tt+tt" decay we define the time-dependent 
B K^W CP-asymmetry 



sin 2a^^ sin Ami, (119) 

where we have used the notation sin2a^^ to make apparent the resemblance with 
Eq. (PPI): we stress however that 2a^^ reduces to 2a when the T^-j^ amplitude 
dominates in Eq. (|39|) , i.e. when the difference between the u- and c- penguins 
dominates over the difference between the t- and c- penguins, which is presumably 



KK 
^CP 



(t) = aff cos Ami - \ l 



„KK 
"dir 
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an extreme case. Conversely, in the absence of long-distance u- and c- penguins, we 
have sin 2a^/ = a^f = f^T], 0. 
Thus, following Eq. (|60D we find 



— COS 2a 



RR 



dir 



cos(2a-2af/) 



(120) 



which reduces to {Pj^oj^l^ = ^k^ro if top-penguin dominates the amplitude. 

Assuming SU(3) and neglecting the (colour-suppressed) electroweak penguin con- 
tributions, we may write 

\P\ = \Pr,oJ^\ . (121) 

ive 

I r - T 9 

= 0. 

'122) 



Thus from Eqs. (^, ( [T^ and QT^ ) we have 
1 - 



a^ir cos(2a — 2a;efr) 



-,RR 
^dir 



RR\ 



cos(2« — 2a^ff' 



Defining the following quantity that can be written in terms of observables 
V = 



a^ir exp{i2acs) - ^^°^° \/l 



-.RR 
^dir 



RR\ 



exp(i 2agg 



= \V\e'^ 
Eq. (|122|) becomes 



cos(2a - ^) 



\V\ 



^RORO 



(123) 



(124) 



As 2aes and 2a^^ are both measured up to a twofold discrete ambiguity, Eq. ( |124| ) 
gives 2a up to an eightfold discrete ambiguity. An explicit example is given on Fig. ||. 



However from the experimental point of view the study of this decay may be 
as difficult as the isospin analysis: First, it is a pure b d penguin decay and is 
thus expected to be very rare (~ 10^^-10^^) [1^. And second, the time-dependence 
of the decay rate may be difficult to reconstruct because the neutral kaons decay 



far away from the primary vertex. This shows the interest of the bound (p7|), very 
symmetrically to the case of the isospin analysis. 



7.4 The Fleischer-Mannel and Marrocchesi-Paver Methods 



Fleischer and Mannel 0, as well as Marrocchesi and Paver |jTT| had already remarked 
that knowing the value of \P/T\ alone leads to the extraction of a. Therefore they 
have used Eq. ([58| ) without explicitly having written it, and without having noticed 
the complete generality of the method. Let us briefly sketch the main points of their 
studies: 
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Figure 9: The eight solutions of the Buras- Fleischer analysis, for the numeri- 
cal example Odir = 0.12, sm2aefr = 0.58, a^/^ = 0.21, sin2a^^ = 0.059 and 
Bj^a-xo/^TT+TT- = 0.058 (see Appendix^, (a) In the {\P/T\,2a) plane, the dots rep- 



resent the central values obtained from Eqs. ( 125 - 12^ ), while the solid curves (limited 
by the small dashes) represent the allowed domain when assuming that Eq. ( \121\ ) is 
affected by a guess 30% relative uncertainty. Some of the solutions merge because 
of this theoretical error, leaving only six separate solutions, (b) The same allowed 
domain is represented in the (p, rj) plane, where it is obtained by plotting the cir- 
cle / jg^ j for 2a varying in the six solution intervals. In the background is shown a 



crude representation of the early-1998 allowed domain pi/. 



Fleischer and Mannel use a first-order expansion in \P/T\. We have shown 
that this approximation, although numerically good, is unnecessary: Eq. 
is exact and not more complicated than its first-order expansion. 



Fleischer and Mannel estimate \P/T\ by assuming Eq. (|109|) and neglecting the 
colour-suppressed contributions to 
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2 




T 







B 



X 



2B^ 



(125) 



while Marrocchesi and Paver use factorization to calculate (in this case \P/T\ is 
just proportional to a ratio of short-distance Wilson coefficients times a CKM 
factor) [|11] 

sin(« + /3)_... (126) 



sin P 



X 0.055. 



The two above equations represent alternatives to the method presented in 
Section ^ although in the second case it is not clear to what extent factorization 
can be used to calculate \P/T\ Note that these two approaches use a single 
model-dependent input, as the method we have proposed in Section (^. 

Both Fleischer and Mannel and Marrocchesi and Paver face the problem of 
knowing \Vtd\ or sin(a + /?)/ sin/5. The first two authors assume simply that 
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|Vtd/(AV^b)| is known from CP-conserving measurements P], while the second 
two authors take the value of [3 as it would be given by the future measure- 
ments of the B J/'^Ks CP-asymmetry and obtain an equation depend- 
ing on a alone |[TT|]. However it is not clear if CP-conserving measurements 
will give |Vtd/(AV^b)| with enough accuracy, and using instead the value of [3 
unfortunately propagates the uncertainty and the discrete ambiguities associ- 
ated with the measurement of (3 into the extraction of a. We have shown 
in § |3.6| that one can avoid these problems by directly writing easy-to-solve 
polynomial equations in the (p, rf) plane, therefore without invoking other in- 
dependent CKM measurements. For the Fleischer-Mannel proposal one should 
write l-P/Tp = A^|l — p — i?7p x Bk-k± / {'^^tt±ito) and report this expression into 
Eq. (|58|) to obtain an equation ^ in the variables (p, without the need to 
know |Vfd/(AVcfe)|. For the Marrocchesi- Paver method one should simply insert 
Rp/Rx = 0.055 in Eq. ( fflD independently of /3. Thus our framework allows to 
improve significantly these proposals. 

• Finally we would like to stress once again the importance of the discrete am- 
biguities. While they are not discussed at all by Fleischer and Mannel 0, we 
believe that the treatment of Marrocchesi and Paver is incomplete: for a given 
value of \P/T\ (inferred from factorization and a given value for (3), they find 
two solutions for a between and vr. We have shown in § |3.4| that there are four 
such solutions which, because of the finiteness of the errors (both theoretical 
and experimental), may merge among themselves. 

8 Conclusion 

We have shown that in the presence of penguin contributions, the information on the 
CKM angle a coming from the measurement of the time- dependent B^{t) ti^tt~ 
CP-asymmetry can be summarized in a set of simple equations, expressing a as a 
multi-valued function of a single theoretically unknown parameter. These equations, 
free of any assumption besides the Standard Model, provide by themselves an exact 
model-independent interpretation of future CP-experiments. 

It is also possible to choose as the unknown a pure QCD quantity, in which case 
the above equations should be expressed directly in the (p, rf) plane, thanks to the 
unitarity of the CKM matrix which predicts relations between the CP- violating angles 
and the CP-conserving sides of the Unitarity Triangle. Whatever the choice of the 
single unknown, as for example the ratio of penguin to tree matrix elements, this 
unavoidable non-perturbative parameter in i? ^ tt^tt' could be compared to Bk in 
the kaon system which allows to report the measurement of ex in the (p, if) plane. 
However the ratio \P/T\ is a much more complicated quantity than Bk, and would 
be very difficult to obtain from QCD fundamental methods. 

^^We have not written this equation, which is not Eq. ([n]), because the Sx^±/(2S^±jro) ratio 
already incorporates a \ V*f,\ factor. 
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Using these analytic expressions, we have assumed some reasonable hypotheses to 
constrain the free parameter. Doing so we have derived several new bounds on the 



penguin-induced shift |2q; — 2Q;efr|, generalizing the result of Grossman and Quinn [15 



One of these bounds is determined by the ratio \'^Bk±-k^ on which one would 

have an experimental value very soon. 

Accepting less conservative assumptions, stronger constraints on a can be ob- 
tained. For example in the limit where the annihilation and electroweak penguin 
diagrams can be neglected, and using SU(3), the knowledge of the — > Kn^ 
branching ratio is a sufficient information to extract the theoretical unknown. As- 
suming a reasonable 30% relative uncertainty (at the amplitude level) on the un- 
avoidable hypotheses, a relatively small allowed domain in the (p, 77) plane can be 
found, independently of any other measurement. This method could be competitive 
with the full Gronau-London isospin analysis, because the latter is plagued by twice 
more discrete ambiguities. From the experimental point of view, our proposal may be 
much more easier to achieve. More generally, if by some other argument a knowledge 
of the modulus of the penguin amplitude — or the ratio of penguin to tree — with a 
~ 30% uncertainty can be achieved, then rather strong constraints on a should be 
obtained. 

However we do not pretend that the theoretical uncertainty on a will be small. 
Rather we believe that this error may be quite well controlled by conservative argu- 
ments. This shows the importance of generalizing our framework to other channels 
sensitive to a: if we are unlucky in the vr+vr" channel, it may happen that we are 
lucky in others. As the problem of the discrete ambiguities is crucial in these analy- 
ses, the modes providing new CP-observables are of a particular interest: for example 
measuring directly the sign of cos 2a;efr rather than determining it from the SM con- 
straints on the UT would be a valuable information, even in the presence of sizeable 
penguin contributions, as it would allow to reduce the discrete ambiguities generated 
when expressing a as a function of the observables and of one model-dependent input. 
It has been shown previously that the analysis of the B ^ pn —>■ 3tt Dalitz plot 
actually leads to the measurement of a kind of cos 2acs (which is of course different 
from the aes in B ^ nn), and we are currently studying the possibility to describe 
this interesting decay similarly to S ^ vrvr Likewise the angular distribution of 
the decay B AA contains also terms proportional to the cosine of an effective a 
angle |8[. 

It is quite clear that all the strategies proposed until now to disentangle the 
penguin pollution in various channels will give different information on a, each relying 
on very different theoretical assumptions and on different observables. Our framework 
allows to treat all these sources of information in a transparent and unified way. Thus 
we will have certainly a strong cross-check of the various results. If this cross-check 
is successful, we may think to combine these results in order to have a more precise 



^^Eventually the _B — > p7r ^ Stt time-dependent Dalitz plot together with the isospin symmetry 
also allows the extraction of penguins |36| . However, such an analysis seems to require a high 
statistics [EtI. 
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knowledge of a. However, we are aware that combining theoretical and experimental 
errors is a difficult problem by itself which is beyond the scope of the present paper. 
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A A Typical Set of Theoretical Parameters 

In this appendix, we define a typical set of parameters in order to compute the relevant 
observables. We assume that 1—6 are exact, and neglect further all annihilation 
diagrams. Thus the amplitudes in Eqs. ( p3l - p5| ), (|39| - PT]) write: 
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(127) 
(128) 

(129) 

(130) 

(131) 

(132) 



Note that in the strict SU(3) limit and neglecting annihilation diagrams, 6' = 6. 



Bt^+t^- (normalization) 


Bt,0^0 


Btt^ttO 




Bk^tt^ 


Bk7T± 


0.75 


0.0455 


0.533 


0.0433 


1.075 


0.948 



Table 1: The various average branching ratios (in units of 10 ^) from our set of 
parameters. They are consistent with CLEO data 



Numerically, we take i3^+,r- = 0.75 x 10~^, which fixes the normalization of the 
amplitudes and choose T real which fixes the origin of phases. Then we choose 
\P/T\ = 0.25 which is a quite sizeable value (cf. § |3.3| ) and 6 = —15° which is a 
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sin 2acs 


"dir 


"dir 


sin2a,f^" 


"dir 


0.117 


0.579 


0.317 


0.209 


0.0592 


0.108 



Table 2: T/ie various CP-asymetries from our set of parameters. 



large violation of naive factorization which gives 6 = 180°. The normalization is then 
given by |T| = 0.826 x v^lO^;B^+^- (in "units of two-body branching ratio"). We 
choose also T^^+t^o = 1.25e~*^° x |T| which takes into account the usual 02 ~ 0.25 
colour-suppression factor and some FSI phases, 6' = +20°, and r„ = 0.3e+*^^° which 
is a ratio of long-distance over short-distance penguin matrix elements. For the CKM 
parameters, we have A = 0.2205 and take p = 0.10, rj = 0.34 which is around the 
center of the early-1998 allowed domain {a = 85.7°) [^. The resulting values for 
the observables are summarized in Tables |l| and |^. Let us stress that these values 
are only indicative and that the real numbers may be very different. Our set of 
parameters results from a compromise between the need to take into account various 
effects in a more or less realistic way and the pedagogical needs (for example it is 
easier to discuss the number of discrete solutions when they are quite well separated, 
which is often not the case in practice). Finally we notice that the penguin- induced 
shift on a is quite large for this set of parameters: 2a — 2aes = +26.7°. 



B Bounds Independent of Direct CP- Violation 

Here our purpose is to derive bounds which are fully independent of a^n, and thus 
are not affected by the experimental uncertainty associated with the measurement of 
direct CP- violation As far as the bound (|16]) is concerned, a different demon- 
stration has been given by Grossman and Quinn ||15|| . 
Consider the bounds (p!6| - |T9D : they all can be written 

= < cos(2a - 2aeff) (133) 

where m is a positive ratio of branching ratios and is expected to be smaller than 2 
(otherwise the bound is useless). If a^ir is not known, then 2aes is not known either. 
Rather one gets from the sin Am t term in (^) the effective angle 



sin 2aes = yl ~ a^^j. sin 2aes ■ (134) 
Since | sin2aefr| < | sin2aefr| one has: 

|cos2Q;efr| < |cos2acfr|- (135) 

As the sign of cos2cIefr is not observable, it can be chosen arbitrarily. It is convenient 
to define 

sign(cos2aeff) = sign(cos2a) , (136) 
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in such a way that ( |135| ) gives 



I cos2aofr cos2a| < | cos2aefr cos 2a| = cos2Q:efr cos 2a . (137) 
Thus Eqs. ( [L33|) and ( |1371 ) imply 



1 — m < Y 1 — cidir '^os ^OefT COS 2a + sin 2aeff sin 2a 

< I cos 2aefr cos 2a | + sin 2aefr sin 2a 

< cos 2aeff cos 2a + sin 2aeff sin 2a = cos(2a — 2aefr) , (138) 
and we obtain the announced result, namely 

< cos(2a — 2aefr) ^ 1 — m < cos(2a — 2aefr) • (139) 



1 '^dir 



It is straightforward to demonstrate an analogous result for the bound (|15|). 
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